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Abstract 

Heavy-light bound states are studied in Lattice QCD with emphasis on parameters of the 
B-system relevant to experiment. Results are obtained on lattices with lattice spacings from 
about 0.15 fm to 0.06 fm corresponding to [3 = 5.74,6.0 and 6.26, and covering sizes from 
about 0.7 fm to 2 fm. From our results at the infinite quark mass limit and from propagating 
heavy quarks with mass of about 1-2.5 GeV we extrapolate to the b-quark mass to obtain 
the decay constant of the B-meson as well as the mass of Aj in the continuum. The neces- 
sary extrapolations introduce rather large errors and lead to the value fs = 180(50) MeV. 
We compare with the results from other lattice groups as well as with predictions coming 
from sum rules. The phenomenological consequences of this value of fs, together with new 
experimental constraints, are briefly discussed. 



1. Introduction 

The extraction of Cabibbo-Kobayashi-Maskawa (CKM) matrix elements from experimental mea- 
surements is hampered by our poor knowledge of hadronic matrix elements. Therefore, in ad- 
dition to increasing the experimental accuracy of the results, we also need non-perturbative 
techniques to evaluate the hadronic contributions. In particular, weak decays of heavy hadrons 
play a key role in probing the hadronic sector of the Standard Model [1] and the big challenge is to 
establish direct CP-violation effects in B decays, as they are predicted in the Cabibbo-Kobayashi- 
Maskawa mechanism [2, 3]. Within this scheme, by virtue of the unitarity of the CKM matrix, 
the CP-violation effects, in principle, can be predicted entirely in terms of CP-conserving exper- 
imental quantities. At present the most important constraints for the CP-violating CKM phase 
S come from the observed CP-violation in the K° — K° system and from results on B° — B° 
mixing and semileptonic B-decays [4, 5]. For the prediction of S, one needs the values of the 
leptonic B-meson decay constant fp and the Bp factor, which are not known experimentally. 
Theoretical predictions for heavy-light bound states like the B-meson, have been notoriously 
difficult and have in the past produced a rather big range of values. They were based on QCD 
sum rules [6], various potential models [7] and lattice calculations at the D-meson mass which 
were then extrapolated to the B- meson, assuming the scaling law [8] 

fp(Mp) \J Mp = constant 
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to hold. Usually a range 100 MeV < fp < 200 MeV had been anticipated in phenomenological 
discussions [3, 5]. 

Needless to say, a full lattice determination of fp, being an ab initio calculation from the QCD 
Lagrangian, would be superior to any other way of predicting this quantity. Practical lattice 
simulations of QCD are faced, however, with a number of problems that must be carefully dealt 
with before one can claim reliable results from this theoretical approach. For B-mesons, due 
to their large mass, the limited lattice resolution (finite lattice constant a ) is potentially the 
most dangerous source of a systematic error: with quarks treated as Wilson fermions, as we 
do throughout this work, such scaling violations are first-order effects in a. In the instance of 
heavy-light quark systems on the lattice, 0(a) terms translate into 0(a Mi) effects and with 
present lattice resolution limited to lattice spacings a > (AGeV) -1 , we have to cope with the 
problem that a Mp > 1. 

In the past couple of years we have learned how to deal in a systematic way with heavy- 
light systems. The infinite heavy quark mass limit is investigated in depth by several groups 
[9, 10, 11, 12]. Smearing techniques have proved essential in extracting reliably the ground state 
with the Fermilab group [13] having developed a method to construct the most appropriate 
mesonic wavefunctions. Allton recently analysed [14] the available static results for the decay 
constant f B tat and found that scaling sets in at (i ~ 6.0. At the same time a number of 
calculations were performed with heavy propagating quarks [11, 12, 15, 16] approaching the b— 
quark from the lower mass side. Investigation of the dependence of fp on the heavy quark mass 
is necessary in order to extrapolate to the B-meson mass. This mass extrapolation, together 
with the continuum extrapolation, introduces the major source of uncertainity on our lattice 
results. Use of an improved action with only 0(a 2 ) effects [17] ablates finite a contamination 
and it is in this direction that future progress is to be expected. 

In this talk we discuss the techniques we used to obtain a lattice estimate for fp as well as for 
the mass of A&. Since this talk is addressed also to non-lattice specialists we start with a brief 
introduction of lattice QCD with emphasis on the relevant features for a numerical calculation. 
We will review the static approximation [18] and compare results from various groups as well 
as with results from QCD sum rules. This will be followed with a discussion of the results 
using propagating heavy quarks and of our procedure to obtain the continuum limit. We will 
also comment on the calculation of the A& mass [19]. In the last section we will see that the 
lattice errors on fp, Bk and Bp are still large and these, in combination with the experimental 
uncertainties, make predictions of the unitarity triangle rather limited. 



2. Introduction to Lattice QCD 

The QCD Langragian density for pure gauge is given by 

c v = -\f; v f^ (2.i) 

where the greek (roman) letters denote Dirac (colour=l,2,...,8) indices and 

F% = %K ~ d,Al + g f cde A^Al 
with f cde the SU(3) structure constants. 
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In lattice QCD we consider a discretized space-time with lattice spacing a as shown in fig. 1. 



n— v 




< — a — > 



Figure 1: Link variables U^(n) defined on a lattice of spacing a 



To put pure gauge QCD on the lattice preserving gauge invariance, the appropriate variables to 
use are the link variables U^(n) [20] 



expiiagoA^in)) 



(2.2) 



which connect side n along the fj, - direction with the neighbouring side n + fi. A plaquette P is 
defined as the product Up of four links round an elementary square. All space is thus covered 
with such plaquettes and the lattice action is the sum of them all. 



Su = P 



E 

p 



1 - ^Tr (U P + U P 



(2.3) 



where 

Up = U^{n)U- v {n + fi)U- tl {n + ft - v)U v {n - v) (2.4) 

That eq.(2.4) reduces to (2.1) can be easily demostrated, by taking the naive continuum limit 
a —7- 0. We have introduced (3 = 6/</q which is an input in a lattice calculation and it is the only 
parameter for pure gauge. It determines the coupling constant and therefore the lattice spacing 
a. The continuum limit is recovered by letting </o 0. Present lattice calculations are done in 
the range 6.0 ~ (i ~ 6.5 which translate into lattice spacing of about (2 GeV ) _1 to (4 GeV) _1 . 



Putting the quarks on the lattice is a much more complicated task that the gauge fields. We will 
not dwell on this problem here but just remark that there is no completely satisfactory way of 
treating fermions, with each known method providing a partial solution to the problem. We use 
the so called Wilson fermions where the doubling problem is avoided at the expence of explicit 
breaking of chiral symmerty [21]. The action with the quarks is then taken to be 



sV 



J2l( x ) M w(x,y)q(y) 



(2.5) 
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= \ q(x)q(x) - k (q(x)(l - !») U^x) c A x + A) + <l{ x + AH 1 + !») Ufa) q(x)j \ 

X I /i J 

with Mw(x,y) the fermionic matrix for Wilson fermions. The hoping parameter k determines 
the naive quark mass via 

m„ = — ( - - — ) . (2.6) 

q 2a \k k c ) v ' 

It is known from chiral perturbation theory [22] that (am^) 2 oc am q and therefore the chiral 
limit is recovered at the critical k value, k c , which is determined dynamically in the lattice 
calculation. In other words breaking of chiral symmetry introduces a fine tuning to be performed 
before lattice results are to match the continuum theory. 

Up to now the theory is defined on an infinite lattice. In numerical calculations one has to 
work on a finite lattice of spatial length L = aNs and temporal length T = aNj with Ns, Nt 
the number of spatial and temporal sites respectively. Therefore in addition to finite-a effects 
one has to consider finite volume contamination. Usually this is checked by performing the 
calculation at the same (3 value with larger Ns- 

The ground state expectation value of some operator O is calculated in the usual way from the 
Euclidean time functional integral 

< O >= lim — / Vq{x)Vq{x)VU ti {x) 0(q,q, U) e -( 5 ^+ 5 f) (2.7) 

T^oo Z J 

which, after intergration of the quarks, yields 

<0>= lim \ [ VU^(x) det[M(U)] 0(M _:l (C/)) e~ Su . (2.8) 

T^oo Z J 

We work in the quenched approximation where internal fermions loops are neglected, i.e. we 
set det[M] = 1 in (2.8). In this approximation the lattice QCD action Sjj couples only next 
neighbours and is thus amenable to Monte-Carlo methods. How good quenching is can only 
be assessed by doing an unquenched calculation and therefore one must bear in mind that all 
quenched results involve an unknown systematic error. At present unquenched calculations 
are starting to yield physical results and one can begin to examine the errors introduced by 
quenching. However the best evidence that the quenched approximation is reasonable comes 
from comparing recent quenched results [23] to experimental values. One finds that low lying 
hadron masses are within 6% ± 8% and decay constants within 12% to 20% of experiment. 



3. Static approximation 



Going from light-light to heavy-light systems the problem we encounter has to do with the finite 
lattice spacing. We would like to study the B-meson and therefore we have to put on the lattice 
a u-quark and a b-(anti)quark with mass of about 4.5 GeV. Since we must keep linear a-terms 
small we require any mass to satisfy 

aM<l (3.9) 

as well as 

1/M < L . (3.10) 
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The second condition is imposed in order that the finite volume effects are small. For the 
B-meson we thus need 



1 



(3.11) 



1/to„ < L 

which is not yet within present computer capabilities. 



Eichten [18] proposed a way to get around this problem. If one looks at the reduced mass, fj,, 
of the B-meson one realizes that it is about equal to the light quark mass to/. Thus fi ~ Aqcd 
which is of course the reason why one cannot apply a non-relativistic treatment. However this 
means that the b-quark mass or in general the heavy quark mass, nih, is much larger that 
typical momenta for the heavy-light bound state. Therefore one can expand in powers of 1 firth- 
Keeping just the zeroth order term in this expansion is referred to as the static approximation. 
In this infinite mass limit the Dirac equation in an external field involves only the the zeroth 
component 

(To-Do - m h ) S (x, y) = 5 4 (x - y) (3.12) 
and solving for the static propagator Sq(x, y) one obtains 



S (x,y) = -i8 3 (x-y) P (*°) 



l + To 



y 0) e -im h (x°-y ) 



(3.13) 



The static (stationary in space) propagator, as expected, has an upper (lower) component which 
propagates forward (backward) in time with an eikonal phase given by the path ordered expo- 
nential 



Pexp 



x 



ig I dz u A c (x, ~ 0uc 

y 



z°)t c 



(3.14) 



This static theory can now be put on the lattice. The heavy quark lattice propagator translates 
into 

S h (x- iy ) = 5 3 (x-y) | e(x°-y )^p-U^x;x°-a)---U^x;y ) (3.15) 



+ 0(y° 



0\ 1 " To 



U (x; x°) ■ ■ ■ U (x; y° - a) 



In eq.(3.15), the exponential prefactor exp(— \x° — y°\ nth) has been dropped, since this corre- 
sponds only to a common shift of all energy levels by the bare mass of the heavy quark. In this 
way the mass of the heavy quark no longer appears and the condition arrih < 1 is relaxed. 

One can envisage going one step further to take into account the first order (1/to/J terms. This 
leads to the familiar non-relativistic Lagragian 

£nrqcd = Q^DqQ - Q 

where Q(x) is a two-component heavy quark; c is one at tree level but gets renormalized, and 
therefore constitutes an additional parameter to be fixed. We remark here that NRQCD is 
non-renormalizable and a finite cut off has to be kept. A number of groups are investigating the 
applicability of (3.16) in the study of heavy-light and heavy-heavy systems [24, 25]. 



D 2 



9o ra 
c a ■ a 



2m h 2m h 



(3.16) 
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Our main approach is to use the static approximation as the leading term in the heavy quark 
expansion and to study propagating heavy quarks with masses less than a -1 i.e. lighter than 
to;,. Interpolating between these results after taking their continuum limit we obtain the real 
estimate for fp [9, 16] or for other quantities of the b-system like the mass of A& [19]. 



4. Pseudoscalar Decay Constant 



The pseudoscalar decay constant fp is extracted from the lattice matrix element through the 

relation 

< 0|A<° o yP >= ^-0Wp/2/ P a 3 / 2 . (4.17) 

Za is the axial current renormalization constant, is the zeroth component of the local 

axial vector current, and \P > denotes the pseudoscalar ground state. The renormalization 
factor Za will be different for propagating heavy quarks and for the static theory. The latter 
we denote by Z s A tat and required a separate perturbative calculation [26]. We define a lattice 
current by 

M ii, & f ) = M*> *) 7o75 l J (x, t) (4.18) 

where 

i J (x, t) = J2 $ J (*, y,W)) Kv, t) (4-19) 

y 

is a smeared light (/) quark field obtained by convoluting with the trial wave function <I> J and 
h(x,t) is the local heavy (h) quark field. 

The aim is to extract the local matrix element by using wave functions optimized to yield early 
ground state dominance in the smeared meson-meson correlator 

(*) = E < M Ls % *) K-y 5 (0, 0)] f > • (4.20) 

X 



A way to extract the local matrix element < OlA^^jP > from the local-smeared, C l °^ , and 
smeared-smeared, , correlator, is to do a fit to the ratio 



R{t)= <0\M 1 ^JP> e~ M ^ 2 . (4.21) 

Ground state dominance is monitored by the plateau in the local mass 



In the static approximation ground state dominance has proven particularly difficult. Unsus- 
pected contamination from excited states was partly responsible for the large values of /^ tat 
obtained at the initial stages. To cure this problem various wave functions have been used by 
different groups [15, 27]. We chose to construct gauge covariant wave functions so that no prob- 
lems due to gauge fixing occur. After experimenting with the parameters and form of these wave 
functions we found that an optimum choice to obtain early plateaus in both smeared-smeared 
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and local-smeared correlators occurs for wave functions with r.m.s. radius of approximately 0.3 
fm. This is a quite reasonable size for a hadronic wave function. In the following we will use the 
'best' wave function, i.e. a 'gaussian type' wave function with parameters n = 100 and a = 4. 
We studied lattices at (3 = 5.74, 6.00 and 6.26. At each (3 value we varied the spatial size to 
investigate the volume dependence. 



4.1 Results in the static approximation 



We consider the quantity 



F = f P ^/M^ 



a s (Mp) 



2/11 



(4.23) 



.a, (M B ) 

which has a well defined limit as Mp — > oo [8]. In fig. 2 we display the dependence of F scaled 
by the string tension a as a function of the spatial length of the lattice also in units of a. 



4.0 




Figure 2: Fa 3 / 4 /^ a * vs La 1 ! 2 at (3 = 5.74, 6.0 and 6.26. The light quark mass is fixed at 
about twice the strange quark mass. The dashed line is the result of fitting all data to the form 
Co - Cexp(-1.5L(7 1 / 2 ). 



From this figure we conclude that for L ~ 1 fm (L^fo ~ 3) the volume effects are small in 
comparison with our statistical errors. Fixing our volume at L ~ 1 fm we then study finite-a 
effects which we expect to be the largest source of systematic error. 



At (3 = 6.0 and using the lattice spacing a =(2.2 GeV) -1 obtained from Allton's analysis [14] 
we display in fig. 3 the dependence of F as we increase the quark mass of the light propagating 
quark from about the strange quark mass to about half the charm mass. 



To obtain the decay constant for the physical mesons we need to extrapolate the light quark 
mass to the chiral limit. This is done in the standard way by first determining the critical k 
value, k c , as explained in section 2. F is extrapolated to k c by fitting linearly to Mp(l,l) for 
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Figure 3: F in GeV 3 / 2 in the static approximation is shown vs the pseudoscalar light-light 
quark mass, Mp(l, /), in GeV 2 at (3 = 6.0, where we took a -1 = 2.2 GeV from ref.[14]. 
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the three lightest quark masses as shown in fig. 3. The results at the chiral limit for our three 
(3 values are shown in fig. 4. 

Following ref. [14] we use for the continuum extrapolation the results at the two largest (3 values. 
This yields F = 0.64(4) GeV 3/2 . If we use the continuum extrapolation of U/y/a = 0.270(12) 
and of F/d 3 / 4 = 1.95(22) (22) from ref. [9, 16] we obtain F = 0.67(7) (7) GeV 3 / 2 consistent with 
that obtained using the a values from Allton's analysis. On the other hand, if one uses the string 
tension to set the scale the value goes down to 0.53(6) (6) GeV 3 / 2 which however is still within 
errors. The value at the B-meson mass translates to 



205 



< ^stat < 

~ Jb ~ 



296 MeV 



(4.24) 



corresponding to the lower and upper values obtained from the string tension and from Allton's 
values of a respectively (only the statistical errors were included in this bound). The Fermilab 
group has recently reported a value of 188(23) (15) (+26) (14) MeV [13], where with their good 
wave functions the isolation of the ground state was optimal. In fig. 5 we summarise the results 
of several groups. Apart from some high values at (3 = 6.0 coming from older calculations, there 
is an overall agreement. In the same figure we display recent results from sum rules [28] also in 
the static limit which show a range of values in agreement with the lattice results. 

The dependence of fg at on the light quark mass is best obtained from ratios where systematic 
errors due to scale and renormalization cancel. Increasing the mass to the strange quark mass 
we obtain in the continuum limit, 
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fstat 
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1.14(11) 



(4.25) 
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Figure 4: F in GeV 3 / 2 in the static approximation at (3 = 5.74, 6.0, 6.26 vs a in GeV 1 taken 
from ref. [14]. The dashed line indicates what the continuum limit is. 
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Figure 5: fs in GeV in the static approximation as a function of the lattice spacing a in GeV 1 
from various groups ( BLS[11], FNAL[13], NRQCD[24], UKQCD[12], APE[10], KENT[29], PSI- 
WUP[9], Sum rules[28]). 
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4.2 Results using propagating heavy quarks 



How good are the results obtained in the static limit? We know that at D-meson mass the 1 firth 
corrections are as large as 40%. In order to answer this question for the B-meson we have to do 
the calculation with propagating heavy quarks on the lower side of the b-quark mass and try to 
interpolate to the static result. For the lattices that we used we can reach masses of about 2.5 
GeV. 

Finite volume effects are on similar ground as with the static theory and will contribute a 4% 
systematic error. To filter out the ground state we will again smear the light quark as we have 
done in the static approximation. Our most serious concern here is the finite lattice resolution. 
When we are close to the continuum limit dimensionless ratios should scale, i.e. should show no 
dependence on the scale. When the dependence on a is weak one may extrapolate linearly as 
was done for the static results. In the fp case however the a-dependence is stronger and such an 
extrapolation is no longer justified. However, what one can still do, is to extrapolate the results 
for lighter masses to the continuum since the a-dependence is weaker there. Therefore we can 
obtain continuum results for pseudoscalar masses in the range 1 — 2.3 GeV. These are shown in 
fig. 6 together with the static continuum result. 
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Figure 6: F in GeV 3 / 2 vs 1/Mp in GeV 1 at the chiral limit and in the continuum. 

We fit the results at finite values of Mp and the static result to a power series in 1/Mp. Using 
c o + C i/Mp + C2/Mp we obtain at the B-meson mass 

f B = 180(50)MeU (4.26) 

This result is in agreement with results from other groups, as shown in fig. 7, where we also 
display results coming from sum rules [28, 30]. 

The dependence on the light quark mass is again best given in the form of the ratio 

f Bs /f B = 1-14(5) (4.27) 

Fig. 8 shows the experimentally measured decay constants for the light mesons and the lattice 
results filling the gap in the heavy-light region. 
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Figure 7: fp in GeV with propagating heavy quarks vs a in GeV 1 . The sum rule results 
were taken from ref. [6,28,30] and are horizontally displayed from a = for clarity (BLS[11], 
NRQCD[24], UKQCD[12], APE[15], PSI-WU[16]). 
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Figure 8: The pseudoscalar decay constant fp in MeV vs the meson mass in GeV. At the Pi- 
meson mass the data referred to as UKQCD[12], BLS[11] and the higher value from PSI-WU[9] 
are static results, whereas the lowest value from PSI-WU[16] is obtained using propagating heavy 
quarks. 
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5. Beautiful Baryons 



We use the same idea of interpolating results from lower mass heavy propagating quarks, that 
we used for fp, to study the baryonic sector. We calculate the mass of A&, a baryon made 
of one b-quark and two light quarks. We avoid computing the mass of the A& directly, but 
rather calculate the mass splittings Aa = M\ b — Mp with respect to the B-meson mass. These 
splittings do not depend on the heavy quark mass in the infinite mass limit and are therefore less 
prone to contamination by finite a effects in the b and c quark mass regions. The interpolating 
field for the A baryon [31] is taken as 

CaW = £((> 6c M*) (u J b (x)C 75 d J c (x))) (e abc h a (0) (^(0)C 75 ^(0)))^ , (5.28) 

X 

where smearing is applied to the light quarks u and d. 

Given the lattice results for the A correlator, C\, and the pseudoscalar correlator, 
C p = J2x < (h(x) 7 5l J (x)) (^0)75/1(0)) >, we perform a direct fit to their ratio 

RA{t) = W) ^ Ae ~ AAt (5 - 29) 

which in the large t limit yields the mass splitting. In fig. 9 we show our results at the three (3 
values together with the result from UKQCD which used the clover action. The a-dependence 
is well within the statistical error and we therefore make a 1/Mp fit extrapolating the low mass 
data to the B-meson mass. On the same figure we display the static results at 1/Mp = 0. The 
highest point is the result of an earlier calculation [31] whereas the lowest is a recent UKQCD 
result using the clover action [32]. The middle point is our result [9]. 




Figure 9: The A& — Mp mass splitting in GeV vs 1/Mp in GeV -1 . The data at (3 = 6.2 come 
from UKQCD [33]. The results at 1/Mp = are static results from [31], [9] and [32] in order of 
decreasing size. 

Using the experimental value Mp = 5.27 GeV we obtain 

M Ab = 5.728 ± 0.144 ± 0.018 GeV (5.30) 
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in good agreement with the result coming from the UA1 collaboration. 



6. Phenomenological consequences 

We briefly comment about possible phenomenological consequences on the expected size [5, 33] 
of CP violation effects in B-decays, using our value of f B = 180(50) MeV. 

As usual, we start from the unitarity triangle 

i=u,c,t 

associated with the b — > d transitions. Using the standard Wolfenstein parametrization we have 
Kt*6 — l^cciKH + Vtd = to a very good approximation. After rescaling by |V C( iV c ^| _1 this defines 
a triangle in the complex plane with vertices (p, tj), (1, 0) and (0, 0) where 

P =(ReV ub )/\V cd V; b \ 1 V =-(lmV ub )/\V cd V; b \ . (6.31) 

We use the updated experimental constraints as given in [34]: 

\V ub /V cb \ = \V us \^p 2 + n 2 = 0.08 ± 0.03 (6.32) 

with |V c6 | = 0.039 ± 0.006 and yj ' p 2 + n 2 = 0.36 ±0.14, 
the K° — K° CP-violating parameter, 

|e| = 4.00 X 10 4 B K \V cb \ 2 \V us \ 2 n [f(x c , x t ) + \V cb \ 2 (1 - p) F{x t )\ = (2.26± 0.02) x 10" 3 (6.33) 

and the B° — B° mixing parameter 

x d = 4.15 X 10" 7 B B T B m- B x f B \V cb \ 2 \V us \ 2 ((1 - pf + r, 2 ) F{x t ) = 0.716 ± 0.04 . (6.34) 

We take t b = 1.63(7) ps and following ref.[34] we assume Bk = 0.8 ±0.2 and take B B = 1. The 
functions F are the QCD corrected versions of the Inami-Lim functions depending on the charm 
mass x c = m^/M^y and the top mass x t = m^/M^y [35]. It is obvious from the expression for 
x d that a large (small) value for f B favours positive (negative) p. 

Fig. 10 shows the allowed (p, r]) range for a top mass nit = 174 GeV, using the above phe- 
nomenological constraints. 

From this figure it is clear that with the present constraints we are not able to pin down the 
unitarity triangle. However there is a good possibility for a positive value of p which will be 
good news for CP asymmetry. 
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m t = 174 (GeV), 130 (MeV) < f B < 230 (MeV) 
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Figure 10: The shaded region is the allowed region for the unitarity triangle. 



7. Conclusions 

Within the quenched approximation we have been able to perform a detailed study of the limiting 
behaviour of fpyJMp by varying the heavy quark mass nih and the lattice spacings a in suitable 
ranges. 

By means of extrapolation, we were able to obtain the decay constant of the B-meson, 
f B = 180(50) MeV, as well as the mass of the A 6 , M Ab = 5.728(144) (18) GeV. Finite lattice 
spacing effects still constitute the largest source of systematic error accounting for most of our 
errors in the lattice results. 

In view of the phenomenological importance of fp to CP violation in the CKM scenario, it 
would be highly desirable to tighten the bounds set in eq. (6.34). This would require pushing (3 
to higher values but in addition using an a-improved action. 
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